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We report on a search for second generation leptoquarks ($2) using a data sample corresponding 
to an integrated luminosity of 110 pb _1 collected at the Collider Detector at Fermilab. We present 
upper limits on the production cross section as a function of $2 mass, assuming that the leptoquarks 
are produced in pairs and decay into a muon and a quark with branching ratio /3. Using a Next-to- 
Leading order QCD calculation, we extract a lower mass limit of M* 2 > 202(160) GeV/c 2 at 95% 
confidence level for scalar leptoquarks with /3=1(0.5). 
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Leptoquarks are hypothetical bosons which carry both 
baryon and lepton quantum numbers and mediate in- 
teractions between quarks and leptons. They appear in 
many extensions to the Standard Model, e.g. GUT, su- 
perstring, horizontal symmetry, compositeness or techni- 
color p] . Leptoquarks which combine quarks and leptons 
of different generations result in flavor changing neutral 
currents, which are known to be highly suppressed [Q. 
While these FCNC constraints do not exclude such lep- 
toquarks, they restrict them to very high masses. For 
example, in the Pati-Salam model j|, the masses are ex- 
pected in the multi-TeV range, and indirect searches for 
such leptoquarks have been made For this search, we 
assume that the leptoquarks couple only to leptons and 
quarks of the same generation. This leads to the classi- 
fication of leptoquarks of three generations, denoted as 
$i, i = 1, 2, 3 in this report. 

In pp collisions, leptoquarks can be pair-produced by 
gluon-gluon fusion or qq annihilation The contri- 
bution to the production rate from direct <i>ql coupling 
is suppressed relative to the dominant QCD mecha- 
nisms M. The coupling strength to gluons is deter- 
mined by the color charges of the particles, and is model- 
independent in the case of scalar leptoquarks. The pro- 
duction of vector leptoquark pairs is also possible. How- 
ever, vector leptoquarks have model-dependent trilincar 
and quadratic couplings to the gluon field 0. In typical 
cases the production cross section is orders of magnitude 



larger than for scalar leptoquarks. The acceptance for 
vector and scalar leptoquark detection is similar, result- 
ing in much more stringent limits on the vector lepto- 
quark mass. 

In this analysis, we report on a direct search for pair 
produced second generation scalar leptoquarks. The pos- 
sible decay channels are: 

^2 ~~ * ''i branching ratio (3 

$2 - ► qiv^i branching ratio 1 — (3 (1) 

where (3 is the branching ratio to charged lepton de- 
cay, and (72 is a second generation quark (c, s). Our 
search for $2 production is based on events having a 
topology including two muons and at least two jets 
(<j> 2 $ 2 — t jj+ y- 2]). The production rate of this decay 
mode is proportional to (3 2 . 

A previous CDF study @ excluded M$ 2 < 
131(96) GeV/c 2 , for /3 = 1.0(0.5) using an integrated 
luminosity of 19 pb~ x . A limit has also been published 
by D0 g, which excludes M$ 2 < 119(89) GeV/c 2 for 
(3 = 1(0.5). Searches at LEP-1 have excluded leptoquarks 
with masses below 45 GeV/c 2 independent of (3 ||. Here 
we present a new limit using an integrated luminosity of 
110 pb" 1 collected during the 1992-93 and 1994-95 Teva- 
tron runs (including the 19 pb _1 of the previous CDF 
study). Searches have also been made for first and third 
generation leptoquark production at the Tevatron [fl0|| , 
LEP § and HERA @. The HI and ZEUS experiments 
at HERA have reported the observation of an excess of 



events at high Q 2 JO] . The interpretation of the excess as 
the production of a first generation leptoquark has been 
ruled out for large (3 by the Tevatron results |l0| . 

The CDF detector is described in full detail else- 
where Jl3[ . Only the detector subsystems that are impor- 
tant in this analysis are mentioned here. The momenta 
of muons are measured in the Central Tracking Cham- 
ber (CTC), a 2.76 m diameter cylindrical drift chamber. 
It is surrounded by a 1.4 T super-conducting solcnoidal 
magnet, covering a pseudo-rapidity (77) range up to 1.1, 
which allows precision measurements of the transverse 
momenta (jpr) of charged particles. Inside the CTC a 
vertex tracking chamber (VTX) allows event vertex re- 
construction using tracks over the range \rj\ < 3.25. Jets 
are detected by the calorimeters, which are divided into a 
central barrel (\rj\ < 1.1), end plugs (1.1 < < 2.4), and 
forward/backward modules (2.4 < |ry| < 4.2). Outside 
the calorimeters, Central Muon drift chambers (CMU) 
in the region |ry| < 0.6 provide muon identification. Out- 
side the CMU lie the Central Muon Upgrade chambers 
(CMP), with additional steel between the CMU and 
CMP detectors to reduce the background from hadrons 
in the muon sample. The region of 0.6 < |?7| < 1.0 is cov- 
ered by the Central Muon Extension (CMX) chambers. 

We use the PYTHIA Monte Carlo generator @ with 
the CTEQ4M parton distribution functions jl6| and 
the renormalization and factorization scales defined as 
Q 2 — Pt' together with the CDF detector simulation 



package, to study the detailed properties of the signal 
for $2 masses between 100 and 240 GeV/c 2 . The signal 
selection criteria are set according to the kinematic dis- 
tributions (e.g. the px of the muons and Et of the jets) 
of decay products determined by Monte Carlo studies, 
optimised to eliminate the background with a minimal 
loss of signal events [[IB) . 

We select events from several different central single- 
muon triggers (ll| with px thresholds of 9 or 12 GeV/c. 
From these events an exclusive dimuon sample is se- 
lected by requiring events with two muons satisfying 
p T > 30 GeV/c (m) and p T > 20 GeV/c (/z 2 ). We do 
not require the two muons to have opposite charge be- 
cause at very high px the charge determination is not re- 
liable. One of the muons is required to have a track from 
the CTC that matches with a stub in the fiducial region 
of the central muon detectors (within 2 cm for CMU, and 
5 cm for CMU/CMP). The muon satisfying this criterion 
is defined as a 'tight' muon. The other muon can be ei- 
ther a tight muon or a 'loose' muon. A loose muon is 
defined as a CTC track that deposits less than 2 GeV of 
electromagnetic energy and less than 6 GeV of hadronic 
energy in the calorimeter tower that it traverses. To en- 
sure good track quality, the track is required to traverse 
at least 75% of the CTC in the radial direction, and be 
matched to an interaction vertex determined by the VTX 
to better than 5.0 cm in the Z direction. Both muons are 
required to be isolated, defined as I < 2 GeV, where I is 
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the sum of transverse energies of all calorimeter towers 
(excluding the one traversed by the muon) within a cone 
of AR — 0.4 around the direction of the muon, where 
AR = -J (A?/) 2 + (A(j)) 2 and <fi is the azimuthal angle. 

The total dimuon identification efficiency, averaged 
over the data sample, lies between 79% at M(<&2) = 
100 GcV/c 2 and 74.5% at M($ 2 ) = 240 GeV/c 2 , with 
the dependence on mass being due to the efficiency of 
the minimum ionizing requirement. The combined aver- 
age identification and trigger efficiency is approximately 
70% over the mass range 100 < M($ 2 ) < 240 GeV/c 2 . 

From this high-py dimuon event sample, we require 
> 2 jets with E^p > 30 GeV and > 15 GeV, respec- 
tively. Jets are reconstructed by an algorithm using a 
fixed cone in rj — (f> space. A detailed description of the al- 
gorithm can be found in Ref . ]l7j . For this analysis a cone 
of 0.7 is used. Both jets are required to be reconstructed 
in the region |ry| < 2.4. Jet energy corrections, due to 
the calorimeter non-linearity, energy deposited outside 
the jet cone, underlying energy from other interactions, 
and the detector geometrical dependence, are applied to 
determine the /i — jet invariant mass. The Z° and other 
resonances such as the J/^ or T are removed by reject- 
ing events with a dimuon invariant mass in the regions 
76 < < 106 GeV/c 2 and M w < 11 GcV/c 2 . After 
applying these requirements, we are left with a sample of 
11 events. 

Cosmic rays can fake high-pr dimuon events; however 



such muons take a finite time to traverse the detector, 
generally entering from the top of the detector and exit- 
ing at the bottom. We use the hadronic calorimeter TDC 
information and a measurement of the opening angle of 
the two muons to reject cosmic ray events. None of the 
11 selected events is identified as a cosmic ray event. 

The numbers of events surviving each selection crite- 
rion are listed in Table |. A total of 11 events pass- 
ing the final selection are shown in Figure [l], plotted in 
the muon-jet invariant mass plane (M^ v.s. Mf^j). From 
two muons and two jets, there are two possible muon-jet 
pairings. We choose the combination having the small- 
est invariant mass difference to determine the leptoquark 
mass for possible candidate events. The reconstructed 
leptoquark candidates of a pair should have equal mass, 
within the experimental mass resolution <j T . 

We therefore search for leptoquark candidates by se- 
lecting events in a 3oy mass resolution region of the 
vs. M^j plane around any given mass, as shown in Fig. ^. 
The mass resolution, estimated from Monte Carlo stud- 
ies, depends on the event geometry and the total event 
energy. Consequently, it varies with the leptoquark mass. 
For example, the maximum values for the mass reso- 
lutions are a r {<P 2 = 120 GeV/c 2 ) = 21.2 GeV/c 2 and 
cr r ($ 2 = 240 GeV/c 2 ) = 46.5 GeV/c 2 . The asymmetric 
mass resolution (oval-shaped regions shown in Figure ^) 
results primarily from the detector resolution, but also in- 
cludes a small probability of misidentifying the jet when 
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additional jets exist in the collision, and cases for which 
a wrong muon-jet pairing combination is chosen. 

The backgrounds for the $2 search include higher or- 
der Drell-Yan processes, heavy flavor decay (from bb or 
tt in the dimuon channel), WW, or Z — > t + t~. An ad- 
ditional background results from W plus multi-jet events 
with a fake muon from energetic hadrons penetrating the 
shielding to reach the muon chambers, or with a hadron 
decay to a muon. These background processes are stud- 
ied using relevant Monte Carlo events samples and ac- 
tual data samples where possible. The major background 
is from Drell-Yan processes (we expect ~ 12 events for 
110 pb _1 ) for which the final state includes a muon 
pair (Z /^y — * plus two or more jets from initial 

or final state radiation. There is a small contribution 
from tt (~ 1.3 events). Other backgrounds are negligi- 
ble due to large muon px and jet Et requirements (bb 
and Z — » t + t~), muon isolation requirements (W plus 
jets), and small cross section (WW). The total estimated 
background is 14 ± 1.8 for an integrated luminosity of 
110 pb^ 1 , before applying the 3 ay mass cut. This mass 
requirement reduces the background substantially, since 
in the background events the reconstructed muon-jet in- 
variant masses are not correlated. For example, we have 
estimated the background contribution to be only 0.3 
events for M($ 2 ) = 200 GeV/c 2 for the 110 pb" 1 data 
sample. 

In the final result, we do not apply a background sub- 



traction procedure, giving the most conservative estimate 
of the cross section limit. The number of expected events, 
N, is given by 

TV = C ■ P 2 • <7(M* 2 ) • Stotal, ( 2 ) 
where C is the total integrated luminosity of the sam- 
ple, P is the decay branching ratio to the charged lepton 
plus quark channel, cr(M$ 2 ) is the cross section for a 
given mass, and Stot is the overall efficiency We evaluate 
the factors entering the overall efficiency as a function 
of $2 mass using actual event samples where possible 
and otherwise simulated event samples. As shown in Ta- 
ble [II], it increases monotonically with M($2); from 9% at 
M($ 2 ) = 100 GeV/c 2 to 22% at M($ 2 ) = 240 GeV/c 2 . 

Possible systematic uncertainties of the measured cross 
section limit have been studied. The major source comes 
from a limited understanding of the initial and final state 
gluon radiation. We have used Monte Carlo samples with 
and without gluon radiation to determine the cross sec- 
tion uncertainty due to this effect. The uncertainty de- 
creases as the $2 mass increases, and it is estimated to 
be 10% for M($ 2 ) = 160 GeV/c 2 . A systematic uncer- 
tainty also results from the Q 2 scale and the structure 
functions used. We compute this effect by varying the Q 2 
scale between 1/4 and 4 of the default value (Q 2 = p 2 -*), 
and by using other structure functions (CTEQ2L fjsf] 
and MRS (A) fisfj ). The jet energy scaling uncertainty, 
which results from detector performance limitations is 
determined by including a 10% energy uncertainty in 



the Monte Carlo reconstruction. Other sources of un- 
certainty, resulting from the detector simulation and the 
limitation of Monte Carlo statistics, are relatively small. 
The uncertainty on the luminosity measurement is 7.2%. 
The total systematic uncertainty varies with $2 mass, 
and is computed to be 15% at M$2 = 120 GeV/c 2 and 
10% at Mj> 2 = 240 GeV/c 2 , as listed in Table || 

We compute the 95% confidence level (C.L.) limits 
on a(pp — > $2^2) 1 , including systematic uncertainties 
with no background subtraction, as a function of the lep- 
toquark mass (see Table [n] and Figure ||). The cross 
section limit for a given $2 mass does not depend on the 
coupling A, and therefore the mass limit does not depend 
on the choice of the theoretical model but only on (3. 
A theoretical Next-to-Leading order (NLO) cross section 
calculation [^0j is also shown in Figure |2|, where the band 
represents the main uncertainty of the calculation com- 
ing from the Q 2 value. Comparing the cross section limit 
to this calculation, a limit of M$ 2 > 202(160) GeV/c 2 for 
13 = 1.0(0.5) is derived. 
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Type of selection 


Number of events remaining 


Total number of sample 


30934 


I s * muon selection (tight cut) 


6844 


2 nd muon selection (loose cut) 


4153 


2 jet cuts 


937 


Jet Et cut 


64 


Invariant Mass cut 


11 


Cosmic ray cut 


11 



TABLE I. The number of events surviving each cut, using 
an integrated luminosity of 110 pb -1 from CDF data. The 
estimated total background from Standard Model sources is 
14 ± 1.8 events. 



$ 2 mass (GeV/c 2 ) 


120 


160 


200 


240 


Total signal detection efficiency,£tot 


0.13 


0.17 


0.20 


0.22 


Systematic error on e to t 


0.019 


0.023 


0.023 


0.023 


Number of candidate events 


1 








1 


Estimated background 


3.8 


1.1 


0.3 


0.1 


a (at 95% C.L.) in pb 


0.34 


0.16 


0.13 


0.19 



TABLE II. Results for different $2 masses for an inte- 
grated luminosity of 110 pb -1 from CDF data. No back- 
ground subtraction is made in the cross section evaluation. 
The candidate satisfying M($2) > 200 GeV/c 2 was previously 
published [6]. 



^350 

-v. 
> 

O 300 



'250 



200 - 



150 - 
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I 1 1 I 1 1 I 1 1 I I 1 1 I I 1 I 1 1 1 I 1 1 I 1 1 1 




/X""*"" 240 




220 \ 




zr — "\ 200 \ 












ffvr\.. 140 1 J. /y 










' , , , , 1 , 


CDF(110pb"') 

, 1 , , , , 1 , , , , 1 , , , , 1 , , , , ' 



50 100 150 200 250 300 350 

MCtij), (GeV/c 2 ) 
FIG. 1. Invariant mass M{p,j) distribution for 
events before applying the mass requirement. A to- 
tal of 11 candidate events are displayed on the M{(xj)2 
versus M(/j,j)i plane. The oval configuration shows the 
limit of the mass requirement for $2 masses between 
100 and 240 GeV/c 2 . M(/ij)i is the invariant mass 
with the higher pr muon, while M(/ij)2 is that with 
the lower pr muon. 
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CDF (110 pb"') 
Theoretical cross section 

CTEQ4 Q = '/ 2 M(4> 2 ) 

CTEQ4 Q = 2 M(* 2 ) 

CDF upper limit 
. 95 %CL 



100 120 



140 160 180 200 220 240 
M(* 2 ) GeV/c 2 

FIG. 2. The 95% C.L. cross section limit for $2 
production for an integrated luminosity of 110 pb -1 . 
The theoretical cross section curves [17] for f3 — 0.5 
and 1.0 are super- imposed. 
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